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Abstract. We present an extensive Quantum Monte Carlo study of the magnetic properties of the mixed- 
spin quantum systems ii2BaNi05 ( R— magnetic rare earth) which show coexistence of 3-dimensional 
magnetic long-range order with 1-dimensional quantum gap excitations. We discuss the validity of the 
performed simulations in the critical region and show the excellent agreement with experimental results. 
We emphasize the importance of quantum fluctuations contained in our study which is absent in previous 
mean-field-like treatments. 



PACS. 7 5.10.Jm, 75.25.-l-z, 75.50.Ee 

I 1 Introduction 

, The competition of different orders in mixed-species sys- 

■ terns is a subject of actual interest, not only in condensed 
I matter physics where for instance mixed-spin systems have 
. been studied intensively in the last years PP but also in 
' other research areas like the physics of ultra-cold atomic 
, and molecular gases where Fermi-Bose mixtures are be- 
ing investigated|2j. The fascination in the mixed-species 

I systems is that they have a very rich phase diagram with 
I coexistence of different orders what leads to interesting 

competing effects. 

Here we want to concentrate on the magnetic behav- 
' ior of a class of mixed- spin quantum systems, R2BaNi05 

(R=magnetic rare earth) which show coexistence of 3- 

■ dimensional (3D) magnetic long-range order with 1-di- 
' mensional (ID) quantum gap excitations [H]. In a previ- 
ous publication we proposed a microscopic model in 
order to describe this coexistence of classical with quan- 

\ turn features in Nd2BaNi05 which we evaluated with the 
1 Quantum Monte Carlo (QMC) method. In the present ar- 
ticle we extend our previous work to the region of critical 
behavior of the whole family R2BaNi05 and present a 
detailed analysis of the QMC simulations which was not 
included in j^. 

The R2BaNi05 systems have two types of spin car- 
riers i.e. S = 1 Ni^+ ions which form antiferromagnetic 
chains running along the a axis of the crystal structure 
and s = i R3+ ions (R= rare-earth), positioned between 
the chains. While perfect isolated S* = 1 chains will behave 
as a Haldane system ^ with no long-range order even at 
r = and a gap in the magnetic excitation spectrum, 
long-range ordering occurs in R2BaNi05 at low temper- 
atures induced by the magnetic s = j ions. Inter- 



estingly, the semi-classical 3D long-range order (i.e. spin- 
waves) coexists with quantum Haldane-gap excitations in 
these systems. This behavior has been intensively stud- 
ied in recent years both experimentally T.'B', as well as 
theoretically 3,4,7 . A first attempt to understand this 
coexistence of classical and quantum orders was done by 
considering a mean field approach PltZ] for the Ni-chains 
interaction which proved to be very meaningful to under- 
stand the basic phenomenon of coexistence. This type of 
approach is nevertheless unable to provide a detailed des- 
cription of the physics involved and one has to invoke the 
use of a microscopic model which should include the main 
interactions responsible for the system behavior. In 0] we 
considered a spin model which describes both the interac- 
tion within the S = 1 chains and the interaction between 
the S = 1 chains and the s = ^ ions in between the chains 
in the following way: 

H = Si^2jSi+l,2j + Jc"^ 5'f 2j(sf 2j-l + Si,2j + l)(l) 

with J > and Jc > 0, denotes spin 1 and s denotes spin 
1/2. The index i runs along the chain direction and j in the 
direction perpendicular to the chains. In Fig. we show 
a schematic representation of the lattice model. We note 
that the coupling between S=l and s=l/2 ions has been 
chosen Ising-like since neutron scattering experiments [H] 
on these compounds show that the excitations associated 
with the rare earths are dispersionless, what indicates that 
the coupling between the Ni and R sublattices must be 
extremely anisotropic and can therefore be approximated 
by an Ising-type term. 

Even very simple, minimal models are hard to evaluate 
and one has to either consider some limiting assumptions 
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Fig. 1. Lattice picture of the model considered for R2BaNi05 
where the big (small) circles correspond to the 5 = 1 Ni spins 
( s = ^ R spins) respectively. J and Jc are the Ni-Ni and the 
Ni-R exchange couplings respectively. 

in order to be able to solve the model analytically or make 
use of numerical simulations, which may prove very accu- 
rate for certain models and properties. In our case, we 
consider the Quantum Monte Carlo (QMC) Method. 

The Hamiltonian has been elaborated in such a 
way that its physical properties are, by construction, in 
qualitative agreement with a wide range of observations 
gathered in the family of rare earth nickelates R2BaNi05. 
It has been shown ^ that quantitative agreement to ex- 
perimental results for Nd2BaNi05 can also be achieved by 
selecting the appropriate value of -y. A first instance of 
this agreement is the ratio between the two basic energy 
scales involved in the problem: i.) The intrinsic gap A of 
the system of independent S—1 chains (when Jc — 0) and 
ii.) the Neel temperature Tpf of the complete system when 
the coupling Ni-Nd is switched on. 

In the well-studied Y2BaNi05 0, the rare-earth 
is non-magnetic and no indication of 3D long-range order 
has been observed, what strongly suggests that the system 
can be described by well isolated S=l Ni^+ chains with 
a Haldane gap A ^ 127K. In the magnetic Nd2BaNi05 
[Tninmn] the Neel ordering is at Tat = 49K . Using 
these quatities we get a ratio r = ^ = 2.59. Analo- 
gously, the gap value of an independent S—1 chain is given 
uniquely in terms of the exchange constant A = 0.410J. 
For Jc = 0.31J, we find that the ordering temperature is 
Tn = 0.163J which gives r = 2.51 (a 3% difference with 
respect to the experimental value). Furthermore, for the 
same value of the transverse coupling Jc — 0.31 J, the stag- 
gered magnetization for the Ni sublattice MNi(r ^ 0) = 
0.79 which corresponds to the value 1.6/zb observed exper- 
imentally. A further validation of this microscopic model 
was obtained from the comparison of the calculated stag- 
gered magnetization as a function of the temperature with 
the experimental one. Actually, the staggered magnetiza- 
tions obtained from the QMC computations for both the 
Nickel (Mni) and rare earth (A/r) sublattices turn out to 



be in very good quantitative agreement with the experi- 
mental results for Nd2BaNi05 OH for Jc = 0.31 J. 

Here, based on the model proposed in our previous 
work we explore the antiferromagnetic/paramagnetic 
transition of the whole family of rare-earth nickelates stud- 
ied so far R2BaNi05 R=Nd,Er,Pr and we propose an ef- 
fective model for the behavior of these materials in the 
critical region. 

We have organized the paper as follows. In section |2] 
we describe the QMC algorithm used for our calculations 
and the data analysis procedure. In section |31 we present 
the critical temperatures of the microscopic model evalu- 
ated with QMC and we compare them with the prediction 
given by the staggered mean field approach proposed by 
Zheludev et al. OE). An effective model for the critical 
behavior is presented in section 0| and the QMC result for 
the spin-spin correlations in the paramagnetic phase in 
sectional 

2 Method and data analysis 

To study numerically the Hamiltonian ^ we have used 
the Loop Algorithm J3] (see also 14 for an extensive 
review), a variant of the QMC method. This method be- 
longs to the family of quantum cluster algorithms and 
provides a very efficient prescription for the sampling pro- 
cess of the configuration space. The key for the success 
of these methods is the following. An accurate simula- 
tion of a spin system requires to gather a large sample of 
(nearly) statistically independent spin configurations in 
which we measure the physical magnitudes of interest. A 
simple procedure to collect a sample of spin configurations 
is by performing local updates involving a few spins in 
each Monte Carlo step. However, if the correlation length 
is sizeable, the number of Monte Carlo steps necessary to 
decorrelate two spin configurations (also called autocor- 
relation time) is large. To avoid such slowing down, it is 
necessary to perform global updates involving clusters of 
spins with a size of the correlation length which: a) pre- 
serve all the symmetries of the hamiltonian, b) keep the 
system near equilibrium. The Loop algorithm gives an ef- 
ficient prescription for constructing such clusters. In this 
way the autocorrelation time is of order one ( i.e. a sin- 
gle Monte Carlo step generates a new spin configuration 
which is nearly independent from the previous one). This 
is especially important since we are interested also in crit- 
ical properties where the correlation length is of the size 
of the system. 

The fact that the coupling between S=l and s=l/2 is 
of Ising type and not Heisenberg simplifies substantially 
the implementation of the algorithm. On the other hand, 
the Hamiltonian ^ is not frustrated and does not show 
the sign problem, therefore all Boltzmann weights appear- 
ing in the evaluation of thermodynamic properties can be 
taken positive after the conventional rotation around the 
z-axis of all the spins in one of the two sublattices of the 
S—1 system. The temperatures of interest for comparison 
with experiment are, as we will see, well inside the scope 
of our QMC method. 
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Fig. 2. Histograms of the magnetization in the ordered phase 
for both sublattices and Jc = 0.31 and T = O.IJ. The black 
(grey) histogram corresponds to the S=l (s=l/2) species. N 
is the number of Monte Carlo measurements of the staggered 
magnetization giving as a result Ms after thermalization. The 
figure suggests that the model has two ground states. 
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Fig. 3. Magnetization in the rare-earth sublattice as a func- 
tion of the inverse size of the system for temperatures very 
close to the Neel temperature Tm ^ 0.183 in units of J. The 
value of the transverse coupling here is Jc = 0.36J. Note that 
the T = 0.184 and T = 0.186 seem to be in the paramagnetic 
phase which is only observable in very large lattices (M_r — > 
as L ^ oo). 



The QMC simulations were performed on finite lattices 
and we carried out a finite size scaling in order to be able 
to compare with experiments and we considered periodic 
boundary conditions Si,+i,2j = ^i,2j and Si^2j" = S^^i. In 
order to illustrate the finite size analysis, in Figure JSJ 
we show the distribution of the staggered value of the 
magnetization after 10^ Monte Carlo steps. 

Since the model has two ground states in the two sub- 
lattices in the ordered phase related by the Z2 symmetry 
(see Fig.(|2Il), the staggered magnetization was computed 
by taking the absolute value of the staggered component 
of the spin operator in each Monte Carlo step in order to 
avoid averaging between configurations around these two 
ground states, and then using the relation 
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where Nmc — 10^ is the number of Monte Carlo steps 
after thermalization. In principle, we considered for this 
extrapolation LxL lattices of size L=8, 16 ,24 , 32 ,48, 64 
spins, where L is the total number of spins including both 
magnetic species. 

In Fig. Q we show the magnetization in the rare- 
earth sublattice as a function of the inverse size of the 
system for a model with Jc — 0.36J where Tjv = 0.183J. 
We observe that for temperatures slightly above Tn, i.e. 
in the paramagnetic phase, very large lattices are neces- 
sary to extrapolate the correct value of the magnetization 
(zero in this case). However, the extrapolation is an issue 
only at temperatures very close to the T/y. As we will see, 
for comparison with experiments we need to compute the 
staggered magnetization at temperatures clearly smaller 
than T/v where extrapolation is straightforward. Actually, 
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Fig. 4. Staggered magnetizations in both sublattices for dif- 
ferent values of L . Far from the critical point, where the cor- 
relation length is small, the extrapolation is straightforward. 
For a given system size each point corresponds to a different 
temperature in the ordered regime T < Tn- 



that is what we see in Figure Q where the relation be- 
tween the staggered magnetizations in both sublattices for 
different L values is shown. In this case the staggered mag- 
netization data may already have converged with lattices 
as small as 32 x 32. 

The Neel temperature was determined by using the 
Binder parameter jl5| g, which is the fourth cumulant of 
the order parameter distribution. 



^3 (M4) 
^ 2 ^ (M2) 



(2) 



At Tpf this cumulant is independent of the size of the sys- 
tem, apart from subdominant corrections to the critical 
scaling. In Fig. ^ we present the typical finite-size scaling 
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Fig. 5. Binder Parameter g as a function of the temperature 
for Jc = 0.36. The hne intersection for different system sizes 
signals the Neel temperature. 

for the Binder parameter for Jc = 0.36J computed in lat- 
tices of sizes L—24, 32, 48, 64 spins where the intersection 
for different system sizes signals the Neel temperature. 
The data for L=8, 16 show consistently significant sub- 
dominant corrections and they were not included in the 
computation of T/v. We suspect that this is a consequence 
of the complex structure of the higher energy excitations, 
probably Haldane excitations in the S=l sector surviv- 
ing in the vincinity of the Neel temperature, what sets 
another lengthscale £,ch ^ 6 slightly smaller than the cor- 
relation length of an independent S=l Heisenberg chain. 
Actually, as we will see below, the energy scale associated 
to these excitations ( the "gap") increases as we increase 
Jc- Therefore only sizes L ^ ^ch should enter the finite 
size analysis [TO] . 

3 Numerical simulations vs. mean field 
approach in the critical regime 

The critical properties deserve special attention because 
the effective staggered field approach proposed by Zhe- 
ludev et al. |SE| -b eing a mean field approximation- can 
break down in the vincinity of the critical region, i.e. at 
the Neel temperature. The idea of the mean field is that 
the behavior of the R2BaNi05 compounds in the ordered 
phase can be described in terms of a 5 = 1 chain of Ni 
ions in a staggered magnetic field induced by the magnetic 
rare earth ions: 

Mm=M{aMK) (3) 

where A4{h) is the staggered magnetization of a S=l spin 
Heisenberg chain as a function of the external staggered 
magnetic field h induced by the s=l/2 ions and a is the 
proportionality constant. Inversely, the otherwise free spins 
s = ^ see the mean-field produced by the neighboring 
S = 1 chains. The staggered magnetization of the R lat- 
tice is then related to the staggered magnetization of the 
Ni lattice by the expression. 

Aff. = Mq tanh(Q;/3MNi) (4) 
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Fig. 6. Functional dependence of the QMC results for Mr vs. 
A^Ni for various values of the coupling constant Jc. The solid 
lines correspond to Eq. 




Fig. 7. Effective coupling a between the two sublattices as a 
function of the Ni-R coupling as obtained from the fits in Fig. 
E|The solid line is a fit to the linear law a = 0.443 Jc + 0.0034. 



where Mq is the effective moment of the rare earth ion 
and (3 = l/keT. 

We aim now to compare our QMC numerical results 
with the predictions for Tat of the staggered mean field ap- 
proach. The coupling between the two sublattices in the 
staggered mean field formalism is encoded in the constant 
a instead of the microscopic coupling constant Jc. There- 
fore the first step is to find a relation between both in the 
antiferromagnetic phase. In Fig. © we present the QMC 
results for Mr — f{Jc(3MNi) for various values of Jc. Note 
that all curves fall on top of each other and, in fact, this 
curve can be described by Eq. where a must be a linear 
function of Jc in the range of transverse couplings studied 
as observed from our QMC results. Furthermore, we have 
extracted the relation a = 0.0034 -I- 0.443 Jc in that range, 
which can be considered linear within our statistical error 
bars (see Figure [Tj). 

At this point it is legitimate to question the relation be- 
tween a microscopic coupling Jc and a frankly phenomeno- 
logical constant a especially at finite temperatures. Since 
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Fig. 8. A/ni vs. the value of the Jc for several values of Jc and 
temperatures in units of J. In the inset Mm vs. the 'effective' 
field on the Ni-subsystem induced by Mr H = 2JcMr |3 for 
various T values. The squares correspond to the DMRG results 
[TT] at T=0. 




Fig. 9. Classical mean-field (solid line) compared to QMC 
results (open circles). Note that a classical mean field theory 
approach fails to correctly relate the staggered magnetization 
in any sublattice. Here we show the best possible fit using the 
rare-earth subsystem. 



our QMC data are collected in the range of temperatures 
Tmm < T < Tn, where T^in was the minimal energy that 
we were able to simulate, it is necessary to study whether 
there is a temperature dependence in a. In the main pic- 
ture of Fig. JSJ we show the relation between Mni and the 
value of the transversal coupling Jc at various tempera- 
tures. Actually the data converge in the T=0 limit to a 
well defined curve MNi(Jc) which at low values of Jc is 
the magnetization curve of the ID S=l Heisenberg model 
at T=0 in a staggered magnetic field computed by Yu et 
al. |17) using the Density Matrix Renormalization Group 
(DMRG) method. 

In the inset of Fig. ^ we directly plot A/ni vs. Hcs 
for various T values assuming that a is temperature in- 
dependent. Now we have that even when the staggered 
magnetization in the R lattice is not at saturation we find 
that the rescaled data are superimposed to the magnetiza- 
tion curve of the ID S=l Heisenberg model at low values 
of Jc- This data collapse is a consequence of the large 
value of the Haldane gap and the low temperatures con- 
sidered. At higher values of Jc, the R sublattice remains 
nearly saturated at the highest values of the T considered 
(T=0.07J), therefore no difference is observed when the 
rescaling is done. The results presented in Fig. (jS)) also 
confirm that the assumption of a temperature indepen- 
dent a is valid in the range of temperatures considered 
and that the mean field approach works very well if the 
chains are weakly coupled to the R magnetic moments but 
there are small deviations at higher values of Jc- Analyz- 
ing the QMC results, the scale that separates the ID from 
the 2D regime is the Haldane gap HoS ~ ^id ~ 0.41J as 
expected. We also note that for the smallest value of Jc, 
Jc — 0.09 J the Neel temperature is so small T/v — 0.06 J 
that some of the points are still in the paramagnetic phase 
MNi(T) = 0. 

In addition, to corroborate the relevance of the quan- 
tum fluctuations in this system, we present in Fig. a 



comparison of a classical mean field model of the staggered 
magnetization in the R sublattice i.e. Mr as a function of 
/3Mr and independent of the magnetization of the Ni sub- 
lattice (solid line) (Brillouin function) with the QMC data 
obtained with the Hamiltonian ^ (open circles). We ob- 
serve that there is a significant deviation between the data 
and the Brillouin function which can be attributed to the 
effect of the quantum fluctuations on the magnetization 
in the R sublattice which are contained in the model Eq. 

but not in a classical mean field model. 

A numerical solution of Eqs. (|3I4|I in the whole AFM 
phase ( < Mr < 0.5 , < M^i < 1.0 ) demands an 
explicit analytic expression for the staggered magnetiza- 
tion as a function of the external staggered magnetic field 
in a S—1 chain. For that purpose we took the DMRG re- 
sults of Ref. ^7] and got the best possible fit. The only 
conditions that we impose to our fitting function are: i) It 
should have a smooth fitting behavior in the wide range 
of values necessary to the numerical solution of H3I4|I ii) 
it should show linear behavior at low external fields and 
iii) it should have asymptotic saturation at high external 
fields. The explicit expression that we obtained was: 

M = Aa.rcta.B.{Bh) + Cta,nh{Dh)+Eh/{l + Eh); (5) 

where A=0. 177457, B=9. 58055, C=0.336168, D=56.5282 
and E=0. 386596. The equations (j3l4|l can be solved now 
numerically confirming in general the excellent agreement 
of QMC and the mean field in the ordered phase but 
there is qualitative disagreement in the T]\[{Jc) relations 
as shown in Fig. H10|l . In Fig. (|10|l we compare the Neel 
temperature as a function of Jc obtained numerically (tri- 
angles) using the Binder parameter as explained in the 
previous section with the mean field solution (dashed line) . 
The most significant feature of the mean field solution is 
its quadratic behavior in a and therefore in Jc. This re- 
sult suggests that the functional relation between a and Jc 
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Fig. 10. Neel temperature as a function of the Ni-R coupling 
Jc in units of J. The triangles are the result of our QMC simu- 
lation in our model Hamiltonian. The dashed line is the mean 
field result given by the equations 14I3II and the solid line is 
the fit to a 2D anisotropic Ising model (see text). 

changes abruptly close to the critical point and a ~ ^fJ^ 
would be more appropiate in that regime. We shall inves- 
tigate this discrepancy in more detail in the next section. 

4 Effective description of the critical behavior 

Our purpose here is to reproduce the relation Tn(Jc) (see 
Fig. HlU|l ) using a model that includes the critical fluctua- 
tions that seem to be missing in the mean field approach. 

To elaborate such a model let us consider the following 
observations: 

i.) The magnetic moments in the rare earth conmute 
with the total Hamiltonian 

K2,±i,H]^0 (6) 

which implies that Mr is a good quantum number of the 
Hamiltonian irrespective of the value of Jc. 

it.) At finite temperatures the model has two ground 
states related by the Z2 symmetry in both sublattices as 
shown in the histograms of the magnetization in the or- 
dered region (see Figure Q where the distribution of the 
order parameter is presented after thermalization.) 

iii.) The spin rotational invariance has been broken 
by the Ising character of Jc- As predicted by Haldane 
|S] and elaborated on microscopical grounds by Gomez- 
Santos the low-energy sector of the XXZ S=l Heisen- 
berg chain can be mapped onto a 2D Ising model. Actu- 
ally, the longitudinal spin-spin correlation function of a 
Haldane chain is identical to the one of the 2D Ising model 
in the paramagnetic phase, including the power-law cor- 
rection to the scaling. The previous relation is an example 
of mapping a quantum model in dimension d to a classi- 
cal model in d-|-z dimensions, where z=l for this case. The 
length of the imaginary time dimension is determined by 
the inverse temperature f3. 

Let us consider then a system of Ising planes coupled 
by the intercalated moments of the rare-earth. In the ma- 



terials of interest we have always T/v < A. In the para- 
magnetic phase Tn < T < A we expect that the autocor- 
relation time in the imaginary direction to be smaller 
than P and the system behaves effectively as a three di- 
mensional Ising model. As the temperature is reduced, a 
crossover takes place when becomes of the order of 
(3. Beyond that point the system belongs to the 2D Ising 
model universality class with different coupling constants 
in the two spatial dimensions. 

The simplest effective model ( i.e. valid in the vincinity 
of the critical point) and compatible with the arguments 
above is an anisotropic 2D Ising model whose critical tem- 
perature is given by the Kramers- Wannier expression |19| : 

sinh(2J|| /Tat) sinh(2 J^/Tw) = 1 (7) 

where Jy = 2.35 and J± = 0.022 are the effective cou- 
plings of the model. In Fig. EH we show the agreement of 
this effective model (solid line) with the QMC data ob- 
tained for Eq. 1^ (triangles). 



5 The paramagnetic phase 

For completion, we present in this section the spatial spin- 
spin correlations in the paramagnetic phase. At T = T/v 
the spin-spin correlation function decays algebraically and 
therefore the correlation length is infinite. As the temper- 
ature increases, the correlation length decreases. If the 
nominal gap A of the independent S=l is large enough, 
the properties of the independent S=l should become vis- 
ible in the paramagnetic phase. The question we are in- 
terested in is whether this crossover is observable in the 
value of the correlation length. 

In Figures and (|12|1 we present the correlation 
functions of the z-component of the spin operators in the 
Ni and R sublattice respectively. We observe that, even- 
though the R-R correlations are two orders of magnitude 
smaller than the Ni-Ni correlations, there are finite corre- 
lations between the s=l/2 spins mediated by the chains. 
We recall that the model doesn't couple the s=l/2 spins 
directly. The inmediate question is: which is the induced 
correlation length between the rare earth magnetic mo- 
ments? 

Following White and Huse [201, we calculated the cor- 
relation length ^ for both sublattices in a 64x64 spins 
lattice by fitting the z-component of the spin-spin corre- 
lation function to the law: 

\{S^Sn\ = AeMj)l-'' (8) 

where A and rj are fitting parameters. In the vicinity but 
not too close to we find a fit of the correlation length 
temperature to a power law in both sublattices. We took 
values in a range ^ ~ 5 — 17 that can be accurately com- 
puted with the procedure described above. The best fit is 
= KiT —Tn)^ where K is the only free parameter^]. 
Such behavior is the one expected for a Ginzburg-Landau 
like description above the critical point confirming that a 



J.V. Alvarez, Roser Valentf: The antiferromagnetic/paramagnetic transition in mixed-spin compounds R2BaNi05 



7 



o.7(r 

0.6 - - 
0.5 © ® 



0.4 - 




Fig. 11. Spin-spin correlation function (z-component) in the 
Ni subsystem and along the chains (Jc — 0.18) in a lattice of 
80x80 sites 




/ 

Fig. 12. Same magnitudes and parameters as in figure IIH 
for the R subsystem. Note the change in scale. 

mean field description is again valid as we depart from 
the critical point from above. Besides, we note that for 
T >> Tn, the correlation length in the Ni subsystem ap- 
proaches the correlation length of a single chain. 

The fitting parameter rj shows a differenciated tem- 
perature dependence in both sublattices i.e 77 = for the 
correlation function Eq. JHJ in the R-sublattice while in 
the Ni sublatticc r] approaches the value 0.5 as the tem- 
perature is reduced in the paramagnetic phase. 



6 Conclusions. 

In conclusion, we have provided a detailed analysis of 
the antiferro-/paramagnetic transition of the R2BaNi05 
mixed-spin quantum antiferromagnets based on QMC si- 
mulations of a lattice model of interacting spin-1 and spin- 
1/2 entities. We give an extensive description of the nu- 
merical calculations and data analysis and comparison of 
our results with experimental observations is very good. 
These calculations go beyond previous mean field type of 



approaches and retain the main features of the behavior of 
the system in the critical region. We propose an effective 
2D anisotropic Ising model to explain the behavior of the 
Neel temperature as a function of the coupling between 
the nickel and the rare earth sublattices. 
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